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I. INTRODUCTION 



Lattice QCD is a first principles approach to calculating low energy hadronic quantities 
using numerical Monte Carlo methods. State-of-the-art calculations are now including 2+1 
flavors of quarks in the weighting of the gauge configurations, thus eliminating the quenched 
approximation. However, partially quenched simulations, where the valence quarks have 
different masses than those of the sea quarks, are still of use when combined with partially 
quenched chiral perturbation theory (PQxPT) Since chiral perturbation theory (%PT) 
comes with a number of unknown low energy constants (LEC's), these LEC's must be 
obtained from non-perturbative methods, e.g., lattice calculations, or from experiment, in 
order to have predictive power. When the number of light sea quarks is equal to three, then 
the LEC's of PQxPT correspond to those of the unitary theory 0, []] , and the LEC's obtained 
from fits to partially quenched lattice data can be used to predict hadronic quantities. Partial 
quenching can therefore be used in order to gain a better handle on chiral fits to numerical 
data, because varying the sea and valence quark masses separately leads to the determination 
of more linearly independent combinations of LEC's. It also allows one to make use of more 
of the available lattice data, since simulating additional valence quark masses is relatively 
cheap compared to generating more ensembles with different sea quark masses. 

In this work we calculate PQ^PT expressions relevant for obtaining K — > irn matrix 
elements from lattice simulations. Although matrix elements of K — > tttt are of importance to 
phenomenology, there are difficulties with extracting multi-hadron decay amplitudes directly 
from the lattice, as expressed by the Maiani- Testa no-go theorem {4]. The implication of this 
no-go theorem is that physical amplitudes can only be computed if the final state pions are 
at rest, or some other unphysical set of kinematics. It was shown by Lellouch and Liischer 
[5I (see also Ref. 6]) that this no-go theorem can be evaded, and that the matrix elements 
can be computed at physical kinematics using finite volume correlation functions. Although 
this method does not require xPT, the physical volume necessary to implement the method 
at physical quark masses is large, and therefore prohibitively expensive given the present 
computational resources. 

An alternative method for calculating K — > tttt from lattice QCD simulations is to obtain 
the leading order LEC's necessary to construct K — > mi from lattice simulations of the 
simpler quantities K — > n and K — > 0. This method was introduced quite some time ago 
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in Ref. [zj]. Given that there are large corrections to kaon matrix elements coming from 
chiral logarithms at higher orders in SU (3) x?T, it is necessary to include next-to-leading 
order (NLO) corrections in the fits to lattice data. This is true both because the light quark 
masses are still relatively heavy in present simulations, and also the physical strange quark 
mass is itself rather heavy. It is an important, and as yet unanswered question whether the 
kaon mass is light enough so that K — > tttt amplitudes can be described by one-loop chiral 
perturbation theory to a useful precision. The issue of convergence is quantity dependent, 
and so must be studied for each quantity of interest. We thus calculate the NLO PQ%PT 
expressions for K — > 7T and K — > matrix elements, including finite-size effects, which are 
needed both to extract LEC's from the lattice, and to assess the convergence of %PT by 
studying fits to lattice data as a function of quark masses. 

In this work we calculate PQxPT K — > 7r and K — > matrix elements in the isospin 
(2+1-flavor) limit. We do not consider the completely non-degenerate quark mass case since 
isospin breaking leads to additional complications [such as (8,l)'s contributing to AI = 3/2 
amplitudes], and these would also not be relevant to current lattice simulations. Thus, 
we restrict ourselves to the 2+1 case in both the sea and valence sectors, but with no 
degeneracies between sea and valence quark masses. We do not present here a complete set 
of formulas necessary to extract all of the NLO LEC's from 2+1 flavor lattice calculations, 
since some of the needed LEC's must be obtained from K — > ixtx amplitudes at unphysical 
kinematics. Even so, the formulas should be useful in extracting leading order LEC's from 
lattice data, and in studying the convergence of the chiral expansion. Note that there are 



many works which discuss the determination of the 



-EC's needec. 
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to construct K — > nn 
HQ, though 
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through NLO in xPT at physical kinematics 0, S Q, Q, Q, Q 

we make no attempt to review the various approaches here. 

For the (8,1) (AI = 1/2) amplitudes there is an additional complication in the partially 
quenched theory coming from the treatment of the gluonic penguin 4-quark operator. For 
the three-flavor theory the situation in PQxPT is simplified significantly if the correspond- 
ing chiral operators are chosen to transform as (8,l)'s under the partially quenched graded 



symmetry group 
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19] . That is the prescription we adopt in the current work. If another 



choice is made, such as, for example, if the chiral operators are chosen to transform under 
the (8,1) chiral symmetry group of the full theory, then additional LEC's enter the calcu- 
lation, making the determination of the desired LEC's more complicated. Although this 
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complication requires some care in the three flavor partially quenched theory, the method is 
still viable, unlike the quenched theory, in which quenched gluonic penguin amplitudes lead 
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19 



20 



22]. 



to large systematic uncertainties 

This paper is organized as follows: in Sec. [IT] we give a review of PQxPT, including the 
effects of the weak Lagrangian, and in Sec. II III we give a quick overview of the calculation 
involved. In Sec. IIVI we review the operator subtraction that is necessary for A J = 1/2 
matrix elements, and introduce the 0( 3 ' 3 ) operator for this purpose. NLO formulas of matrix 
elements of this operator are calculated for use in later sections. We present results for the 
(8,8) electroweak penguin operators for the K —> and K — > n processes in Sec. |V] where 
we also give the physical K — > -kit amplitudes for completeness. In Sec. |VTJwe present results 
for the (27,1), AI = 3/2, K — > it matrix element, and in Sec. IVHI we present the results for 
(8,1) and (8,1) + (27,1) operators for K — > and K — > it, including the operator subtraction. 
In Sec. IVIIII we discuss the finite volume corrections for the results presented in this work. 
We conclude in Sec. IIXI and include relevant function definitions and the chiral logarithm 
contributions in a set of appendices. Appendix [H] provides an erratum for Refs. [10] and 

Hi- 



ll. PARTIALLY QUENCHED CHIRAL PERTURBATION THEORY 

We use the standard formulation of partially quenched chiral perturbation theory 



(PQxPT) introduced in |23|, |24|. In this formulation, the valence quark loops are removed 
by introducing "ghost" quarks with the same masses and quantum numbers as their va- 
lence counterparts, but which obey opposite statistics. The chiral symmetry group for a 
partially quenched theory is graded; in general one takes it to be SU (N va i + iV sea | N va {) L (g) 
SU (N va \ + iV Bea | N va _i) R . For the purposes of this work, we set N va x = N sca _ = 3. Specifically, 
we have three valence quarks denoted as x, y, and z\ three sea quarks denoted as d, and 
s; and finally three ghosts: x, y, and z. 
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A. Strong Lagrangian in PQxPT 



As explained in Ref . 11[ , in the partially quenched theory, operators are written in terms 
of the chiral field 



exp 



2z$ 
T 



(1) 



where / is the meson decay constant in the £77(3) chiral limit (normalized such that the 
physical f n « 130.7 MeV), $ is a 9 x 9 matrix containing the meson fields, 




(2) 



where is a 6 x 6 matrix of pseudoscalar mesons constructed out of valence and sea quarks, 
<fi is a 3 x 3 matrix containing mesons constructed with two ghost quarks, x (x ) is a 3 x 6 
(6 x 3) matrix containing fermionic mesons made out of one quark and one ghost quark. S 
transforms under the graded chiral symmetry group as 



(3) 



with L G SU (6 | 3) L , R G SU (6 | 3) R . Operators in the chiral effective theory are con- 
structed from the quark-level operators out of X and other objects (such as the quark 
charge matrix and mass matrix, for example) such that they transform the same way under 
the chiral symmetry group. 



The leading-order (LO) strong Lagrangian is given by 25] 



C% = J -str [d,E^£t] 
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str [MY* + EtA/ft] 



(4) 



where the superscript (2) indicates that this Lagrangian is valid to O (p 2 ) in the chiral power 
counting scheme, and M. is the quark mass matrix 



M = diag (m x , m y , m z , m u , m d , m s , m x , m y , m z ) 



(5) 



Note that this corresponds to the quark vector composed of valence quarks, sea quarks, and 
ghost quarks 

q = (x,y,z,u,d,s,x,y,z) T . (6) 
The supertrace is defined as follows: for a 9 x 9 matrix 

^6x6 -^6x3 
^3x6 -^3x3 



U, 



9x9 



(7) 



in which sub-matrix A is the top-left 6x6 diagonal block and D is the bottom-right 3x3 
diagonal block, then 

str (U) = tr (A) - tr (D) . (8) 

We set the valence x and y quark masses equal, and we set the sea u and d quark masses 
equal, 

m x = m y , m u = m d . (9) 

Thus we work in the isospin limit in both the valence and sea sector, and we present results 
for both this (2+1-flavor) case and the 3- flavor case (degenerate valence quarks). 

At NLO in the full theory [O (p 4 Y|, t he strong Lagrangian involves 12 additional opera- 



tors with undetermined coefficients 



25 



26 L There is an additional 0(p 4 ) operator which 



appears in the partially quenched theory [27] . though this operator does not contribute to 



the quantities considered in this work. The NLO operators of the strong Lagrangian relevant 
for the current work are 

( 4 st) = str [L 2 ] str [S] , 
{ 5 st) = str [L 2 S] , 
Of ) = str [S] 2 , 

Oi st) = ^str [S 2 - P 2 ] , (10) 

where 

S = 2B (A^S 1 " + £M) , 
P = 2B (A^S 1 " - ZM) , 

L^iZd^. (11) 

As follows from the strong Lagrangian above, the leading-order mass of a bare pseudo- 
scalar meson is 

m|- = B (m; + mj) , (12) 

where is the mass of meson rrii and rrij are the masses of the quarks qi and qj 
can refer to the sea, valence, or ghost quarks in this case). In our partially quenched 
amplitudes we assume that the light quark masses are all light enough compared to the rj' 
mass so that the r/' can be integrated out. As demonstrated in Ref. [sj, this is the case where 
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the LEC's of the partially quenched theory with three sea quarks correspond to the LEC's 
of the unitary theory. 

In the following we adopt the notation that the masses of mesons which are constructed 
out of two different flavors of quarks are labelled in terms of their quark consituents, re- 
gardless of whether they are sea or valence, for example m xy or m zs . For any flavor- neutral 
meson, we use = rridd or m>x = m xx for mesons in the "flavor basis." Due to the 
disconnected propagators which arise in the flavor-neutral sector, this is distinct from the 
"physical basis," where the relevant mesons are the 7r° and 77. These only arise in the sea 
sector, and we will use the fact that 

m lo = ™ 2 d , 
m l = \ ( 2m s + m V) , 

in the isospin limit (m u = rrid)- 

The propagators for flavor neutral mesons are obtained by following the prescription in 
Ref. [s| in Minkowski space: 

jT^ = ggfj i (P 2 - ™d) (p 2 ~ m s) / 13 x 

" p 2 - m| + ie 3 (p 2 - m|) (p 2 - m 2 j3 ) (p 2 - m 2 ) 

The propagators for flavor off-diagonal mesons are 

= 2 l % T ■ (14) 



where 



1 j E {x,y,z,u,d,s} 
e 3 ={ . (15) 

-1 j e {x, y, z} 



B. Leading-Order Weak Lagrangian 

In full xPT, we group the weak operators appearing in the K tttt transition by their 
chiral transformation properties in the SU (3) L ® SU (3) R symmetry group. We can carry 
this same idea over to PQxPT, where we extend their definition into the graded group using 
the expanded chiral field E and replacing traces with supertraces. Except for those cases 
discussed explicitly (such as the quark charge and mass matrices, for example), when going 
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from unquenched to partially quenched x?T, operators are replaced as follows: 

(16) 



^3x3 03x( 
06x3 ®6xt 



where the upper left block of this matrix is the 3x3 block corrsponding to the valence 
sector. 



The leading order weak operators are 7|, |8|, |9|, [ll|, [28] 



Ofo 8) = str [A 6 £Q£t] 
o££i = str [Ae^E^Et] 

= 4 (S^Et)J (E^Et); (17) 

where Q is the quark charge matrix, (Xe)ij = ^a^j2i and the tensor t 1 ^ is symmetric on any 
indices and traceless on pairs of upper and lower indices, and its elements are chosen to 
pick out the AS = 1 transitions; it thus plays a similar role to that of A6 for the other 
operators. However, we will defer the actual determination of its non-zero elements until 
Section IVI Al where we will use this tensor to further divide the operator into the isospin 
3/2 part and the isospin 1/2 part and directly evaluate their respective amplitudes. The 
isospin decomposition of K — > ir matrix elements is given in Appendix lAl 

There is a choice to be made for the quark charge matrix, Q, above, which enters in the 



electroweak penguin operators [ll| . We could either assign charges to ghosts such that they 
cancel out the electroweak valence quark loops, or we could make them uncharged. In this 
paper we derive the amplitudes with the electroweak penguin operators for both choices, 
which we denote as Qi and Qi- We always assign zero charge to the sea quarks, since that 
is what is typically done when generating lattice gauge fields. The two choices of charge 
matrix are: 

g 1 = diag(2,-l,-l,0,0,0,2,-l,-l), 

Q 2 = diag (2, -1,-1, 0, 0, 0, 0, 0, 0) . (18) 

We discuss the weak operators which contribute to K — * it and K — > tttt at next-to- 
leading order in subsequent sections. 
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(a) (b) 

Figure 1: Diagrams contributing to K — > at NLO. The gray square is the insertion of a NLO 
weak vertex, and the small dot is an insertion of the LO weak vertex. 

III. DETAILS OF THE CALCULATION 

To make complete use of lattice data in extracting LEC's relevant for K — > hit, it is 
important to work in the non-degenerate m x = m y ^ m z case. Since the K — > ir amplitudes 
do not conserve 4- momentum for m z ^ m x , the weak operator must transfer a 4 momentum 
q = p xz — px- In our calculations we restrict ourselves to the case where both initial and 
final mesons are at rest, so q = {m xz — mx, 0, 0, 0). 

The NLO diagrams contributing to K — > and K — > it are given in Figs. [1] and [2, 
respectively. The external legs are always mesons made of two valence quarks, while the 
internal loops in the partially-quenched theory consist of valence-ghost, valence-sea, and 
valence-valence mesons. In addition to these diagrams, the renormalization of the external 
legs (wave-function renormalization) via the strong interactions must be taken into account. 

The logarithmic expressions presented in the appendices of this work are quite lengthy. 
Thus, checks are necessary. The first check was that the one-loop insertions cancel those of 
the divergent counterterms, and this check was performed for all expressions in this paper. 
Another check is that an expression reduces to some other in the appropriate limit. All of 
the logarithmic expressions in this paper reduce to those in Refs. in the ap- 

propriate degenerate sea quark and full QCD limits. 1 Finally, all one-loop expressions in this 
work were computed separately by at least two of the authors, using independently written 



code. The codes used were the FEYIN 



CALC package [29] written for the Mathematica [30] 



system, and the FORMCalc package 31[, which interfaces FORM 32] with Mathematica 



1 Note that Ref. [8| contains errors that are corrected in Ref. [9( , and we agree with the latter 
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(a) (b) (c) 



Figure 2: Diagrams contributing to K — > it at NLO. As in Fig. [TJ the gray square is an insertion 
of a NLO weak vertex, and the small dot is an insertion of the LO weak vertex. The large dot is 
the insertion of an 0(p 2 ) strong vertex. 



IV. SUBTRACTION OF AI = 1/2 AMPLITUDES 



In general the AI = 1/2 matrix elements of four-quark operators have a power divergent 
part due to the four-quark operators mixing under renormalization with lower dimensional 
operators when using a lattice regularization. This power divergence reduces to a quark 



28|. Following Refs. 



in order 



bilinear times a momentum independent coefficient 
to remove the power-divergence of A J = 1/2 operators, we perform a subtraction using the 
dimension three quark-level operator, 



e (3,S) sj-^.^J 



(19) 



This subtraction must also be performec 



n 



in PQxPT for comparison with the subtracted 



lattice results. Again following Refs. [7], 
corresponding to the (3, 3) operator in Eq. ( fl~9]) is 

egf =«( 3 < 3 >str [A 6 S] 



28| . the lowest order [(9(p )] chiral operator 



(20) 



where the low-energy constant c^ 3 ' 3 ** can be related to the coefficient of the mass term in 
the leading-order strong Lagrangian 111 ]. 

f 2 B 



(3,3) 



As explained in Ref 28|], the 0( 3 ' 3 ) operator can be used to remove the power divergences 
to all orders in the lattice calculation. This subtraction is performed to NLO in PQ^PT 
explicitly in the sections that follow. To this end, we require the higher order chiral operators 
of G>( 3,3 ). The terms up to 0(p 2 ) needed for this work are 



0(3,3) = @ 



(3,3) 
LO 



I> 3 ^ 



(21) 
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where i takes the values 4, 5, 6, 8, H 2 , and where 

1 



0' 4 = -str [A 6 £] str [d^d^] , 

C 6 = 2fi str [A 6 £] str [M f £ + £ f M] 

Cg = 2fi str [AeEA^E] , 

C^ 2 = fiostr [X 6 M] . 



(22) 



The coefficients c 33ii of the operators 0[ are related to the Gasser-Leutwyler coefficients by 
c 33,i = —&B Li, a relation similar to that for the leading order coefficient, a^ 3 ' 3 \ given above. 
To NLO, the K — >• matrix element for 2+1 valence flavors is 







0(3,3) 



2% 
1 



a 



(3,3) 



1 



+ 



-5Z T . 



1 + R x (m v , m z ) - R v (m x , m x ) 



+ 



1 + Rz(m v , m x ) - Rr,(m z , m z ) 



£(m 2 z ) 



- U{m 2 xd ) - U{m 2 xs ) - U{m 2 zd ) - U{m 2 zs ) 



+ 



Rr,{m x , m x ) + R v (m x ,m z ) + R v (m z , m z ) 



- R x (m v )£(m x ) - R z (m v )£(m 2 z ) 
~ ^y^{ L 8™L + L 6 (2m 2 D + m|) 



(23) 



where the chiral logarithms £{m 2 ) and £(m 2 ) are defined in Appendix [BJ along with the 
residues R x (m a ), i?a.(m a ,mb). The wave-function renormalization 5Z XZ is given in Ap- 
pendix o 

To NLO, the K — > ir matrix element (also for 2+1 valence flavors) is 



7T 



0(3,3) 



K~ 



2 ,„ s , I6B0 



f 2 



L 5 m x m xz - 2L 8 [m 2 x + m 2 xz ) 



- 2L 6 (2m 2 D + m|) \ + (vr 



0(3,3) 



logs 



(24) 



For clarity, the rather lengthy logarithmic contribution is given in Appendix [Dl For degen- 
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erate valence masses (m x = m y = m z ), the K — > tt amplitude simplifies to 

2 



7T 



0(3,3) 



\ dcg.val. 



(3,3) 



[1 + sz x ] 



+ 



m 



x 



R x (m v )£(m 2 x ) - 2R n (m x , m x )P(0, ml, m 



,,u x ) 



+ 



-1 + R v (m x ,m x ) 



£{m\) + 2£{m 2 xd ) + £(m 



- R ri (m x ,m x )l(m 2 n ) 



+ 



2m 2 x [ 1 - R v (m x , m x ) ) + R x (m v ) 
P 



i(m 2 x ] 



(L 5 - 4L 8 ) - 2L 6 (2m^ + m|) |, 



(25) 



where £(m 2 ) and /3(g 2 , m\, m 2 ) are defined in Appendix iBl and 5Z X is given in Appendix O 
These expressions are used below when performing the power divergent operator subtractions 
that are necessary in order to obtain the physical amplitudes in which we are interested. 



V. WEAK MATRIX ELEMENTS WITH (8,8), AI = 3/2 AND 1/2 OPERATORS 

In this section we present the results for the chiral operators which transform as (8,8) 's 
under the chiral symmetry. These correspond to the electroweak penguin 4-quark operators. 
Formulas are presented for K — > and K — > ir for nondegenerate (m x = m y ^ m z ) valence 
quark masses, as well as for K — > n with degenerate valence quark masses. The power 
divergent subtraction is discussed for the AI = 1/2, K — > tt amplitude. Since K — > and 
K — > tt are sufficient to construct K — > txtx to NLO at physical kinematics for the (8,8) 's, 
we present the physical K — > tttt amplitudes as well. 
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Following Ref. 11], the form of the operator CK 8,8 ) through NLO in PQyPT is 



(9 (8 ' 8) =a 88 str [AeSgst] 

+ c 88 , lS tr [XeL^Q^L^ 
+ c 88 , 2 str [A 6 L M ] str \LQ0lf\ 
+ c 88 , 3 str [\ 6 {ZQ^,L 2 ]] 
+ c 88i4 str [A 6 £}] 
+ c 88i5 str [A 6 [£Q£ f ,P]] 

+ c 88 , 6 str [A 6 Sgst] str [S] , (26) 



with S, P, and as defined in Eq. (fTTj) . 

A. isT — > amplitudes for 2+1 valence flavors 

As explained in Section III there are two choices for the quark charge matrix for the 
operators in the (8,8) representation. If we set Q = Qi, we obtain for the K — > amplitude 

(0 |0< 8 ' s '| K") Qi = ^ass [-21 (my - I (ml) + 21 (ml,) + I (ml)] 

K - 4) ■ (27) 



4 



/ 

If we set Q = Q2, we obtain 

<0 K°) Q2 = j 3 a S8 {-£ (m\) + 2i (ml) - i (m|) - 2^ «) " « «) 

+ 2^ K) + « (™L) } - yc 88 , 4 (m^ - m y . (28) 

B. K — > 7T amplitudes for 2+1 valence flavors 

The process K — > 7r must be separated into its A J = 3/2 and A/ = 1/2 pieces, and we 
;ive the explicit isospin decomposition in Appendix |A] For the AI = 3/2 amplitudes we 
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have 



(tt + \O^W\K + ) Q2 = (7r + \0^ 3 / 2) \K + ) Qi 

4 



p +p 



1 + c 88 , 2 ) m xz m x 



+ 2 (c 88j4 + c 88j5 ) (m 2 xz + m x ) + 2c 88j6 (2m? D + m|) 

+ ( 7r +IO(8.8)(3/2)| jK -+\ 

\ I I / Qi.logs 



(29) 



For brevity, we give only the analytic part of these matrix elements here; the logarithmic 
contributions are given in Appendix [El 

For the AI =1/2 amplitudes, we are ultimately interested in the subtracted versions, as 
discussed in Sec. IIVI We expand the amplitude (0 |G( 3 > 3 )| K ^ to leading nontrivial order, 
and take the ratio 



[0\O^\K°) Ac 8 



-B (m z - m x ) + 



(logs) 



(30) 



(0\Q( 3 ^\K°) 0,(3,3) " uv "" 2 -" x > ' /2 a (3,3) 

where higher order terms in chiral perturbation theory are omitted. The power divergent 
contribution is proportional to m z — m x , and this is true to all orders in the chiral expan- 
sion by CPS symmetry 3^|. Thus, the ratio of LEC's containing the power divergence, 
— 4c 88 4 i?o/(a^ 3 ' 3 ' ) ), can be extracted from the corresponding lattice matrix elements, since 
the mass dependence of the divergent piece is known to all orders of the chiral expansion. 
We perform the operator subtraction using this ratio and the amplitude (V + |0( 3 > 3 )| K + }, 



7T 



o 



(8,8)(l/2) 
sub 



<7r + |0 (8 < 8)(1/2) |^ + ) 
, 4c 88 , 4j B (m 



Q 



a 



(3,3) 



7T 



6 



(3,3) 



K~ 



(31) 



where by CPS symmetry the power divergence is removed to all orders in x?T. Through 
NLO in xPT we have, 



7T 



o. 



(8,8)(l/2) 
sub 



K 



,7T 



| (8,8)(l/ 2 )|^ +)Q _^ C88i4mL 



(32) 



These relations hold for either Q = Qi,Q2, and lead to 



7T 



o. 



(8,8)(l/2) 
sub 



K~ 



P 



j(-c 8 8,i + c 88 ,2 + 2c 88i3 ) m xz m x + 4c 88j4 (m 2 xz + m 2 x ) 



+ 4c 88i5 (m 2 xz + m 2 x ) + 4c 88>6 (2m 2 D + m|) 

+ / + | (8,8)(l/2)| K+ \ 

\ i i /<yi,iogs 



(33) 
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7T 



o. 



(8,8)(l/2) 
sub 



Ti 



o, 



(8,8)(l/2) 
sub 



+ -^<!3- 



m 7 



t{m 2 x ) +£(m%) -2£(m 



2 ) 

xz J 



+ 6m xz m x (P(q 2 , m 2 xz , m\) - p(q 2 , m 2 xz , m|)) 



(34) 



The logarithms appearing in Eq. ( l33l) are given in Appendix [El 



When we have degenerate valence quark masses (m x = m y = m z ), the above formula 
can be simplified. However, for some terms, especially those which involve residue functions 
Rx { m ai m b)i taking this limit is non-trivial. Thus, we give the degenerate valence K — > n 
amplitudes explicitly. (The degenerate valence K ^ matrix elements vanish due to CPS 



symmetry 



34|.) The subtracted amplitude for the degenerate case is given by 



7T 



o. 



(8,8)(l/2) 
sub 



K +\ dCg ' Val - = / Tr +| g (8,8)(l/2)|^+\dog-val. | 4c 88 , 4 (2£? m a ,) , 



0(3,3) 



K~ 



(n+\0^ 1 ^\K + ) 



deg.val. 



-pC88A m X ) 



(35) 



where again the second equality is correct through NLO in x?T. In the degenerate case the 
amplitudes are the same for Qi and Q%, 



/ n + | (8,8)(3/2) I K+ \ d ^ = 1^8 (1 + SZx 



P 
+ 



f 4 



[21 (m 2 xd ) + I (m 2 xs )] + 2m 2 J (m 2 x ) 



+ -p [(-Cb8,i - c 8 8,2 + 4c 88i4 + 4c 88i5 ) m\ 



71 



o 



(8,8)(l/2) 
sub 



deg.val. 8« 88 



Ql,2 



+ 2c 88j6 (2m 2 D + m 2 s )]. 
(1+6Z X ) 



(36) 



/ 
+ 

+ 



f 4 
4 



-|[2£(my+i(mL)]-m^K) 



~ c 88,l + c 88,2 + 2c 88j3 + 8c 88j4 + 8c 88]5 ) 

+ 4c 88j6 (2m|, + mg) 



(37) 



C. K — ► 7T7T amplitudes in full QCD 

The LEC's needed to construct the (8,8), K — > irn amplitudes at physical kinematics 
through NLO can be obtained from the K — > ir and K — > amplitudes given above. The 
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extraction of LEC's is essentially unchanged from the case of 3 degenerate sea quarks treated 



in Ref. 



For completeness, we present the physical K — > titt am plit udes, which were 



calculated originally in Refs. [8|, [9J , and subsequently checked in Refs. 

4za 88 Yli 



11 



35|. 



(n + n-\O^ 3 ^\K ) 



+ 7 — H [(~ c 88,2 ~ c 8 8,3 — 2c 88]4 — 2c 8 8,5 — 4c 88)6 )m 2 ^ 



(-C88,i - c 88 , 2 + 4c 88j4 + 4c 88j5 + 2c 88i6 K] + (7r + 7r-|0 (8 ' 8) ' (3/2) |A°) log s 



(38) 



(tt+tt-I^ 8 ' 8 )^ 1 / 2 )^ ) 



Yli 



[(-Cs8,i - c 88,2 + 4c 88j4 + 4c 88j5 + 8c 88i6 )m 



+ (-c 88 ,i + c 88j2 + 2c 88j3 + 8c 88i4 + 8c 88i5 + 4c 88j6 )m 2 ] + (n + n~\O i8 ' 8Ul/2) \K°) 



K 



logs- 



(39) 

The logarithmic terms are given in Appendix [El Note that in Eqs. (|38|) and (139^) . the decay 
constants appearing in the tree-level terms are the physical decay constants (correct to one 
loop). When constructing K — > titt amplitudes using Eqs. (|38l) and (1391) . one should use the 
physical decay constants in the tree-level expression, as determined from lattice calculations 
or experiment, in order to avoid double counting a subset of the one-loop corrections. 



VI. K ^vr FOR THE (27, 1), AI = 3/2 CASE 

The operators which transform as (27,1) 's under the irreducible representation of the 
chiral symmetry group give the dominant contribution to Ke(A2), i.e., the real part of the 
AI = 3/2, K — > tttt amplitude. In this section we review the (27,1) chiral operators that 
are needed through NLO, and we give results for the NLO AI = 3/2, K — > ir amplitude. 
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A. Definition of the C( 271 ) operators 

Following 0, [ill, ^ ne operator in the (27,1) representation can be written as 

+ c^lg (S)f 

+ c 27>2 T% (P)f (P)J 

+ cs»r..2S (5)? 

+ 0*7,24?$ (WW)f (Wnj (40) 

where S", P, L M are defined in Eq. (fTT!) and 

= 2 (9^ + d v Lfj) . (41) 

The tensor has different elements depending on which isospin we are projecting. To 
project the AI = 3/2 operator, we set 

rpl3 rpSl rplS rpSl 

1 12 ~ 1 12 — 1 2l — 1 21 — 2 

m = m = -~ , (42) 

whereas for the AI = 1/2 operator, we set 

rpl3 rpSl rpl3 rp'Sl 

1 12 ~ 1 12 — 1 2l — 1 21 — 2 
J 22 — J 22 — 1 

TM = 7f = ~\ ■ (43) 

In order to adapt Eq. (140j) to the partially quenched theory, we must promote T to a 9 4 
element tensor, although many components will remain zero (only the 3 4 block corresponding 
to the valence quark sector will have non-zero elements). To take into account the graded 
nature of the group, we multiply by factors of £j defined in Eq. (Tl~5|) . such that 

o^=x;^ii]g(E^ t )f( s ^ Et ); > (44) 

ijkl 
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where we display the summation over i,j, k, I explicitly for clarity. 

There is another equivalent approach to obtaining the partially quenched operators for 

to rewrite Eq. (1401) in 



the (27,1) case. It is possible, as illustrated in Appendix D of 
terms of traces over the various operators. The partially quenched theory is then obtained 
in the usual way by changing traces to supertraces, and we obtain 

(27,i),(3/2) = s t r [A 6 £0 /i £ t ]str[A£0"£ t ] +str[A 3 E9 M S t ]str[A 4 i:^S t ], (45) 

(27,i),(i/2) = stT [x 6 Ed^]str[BEd^] + str[A 3 £<9 /1 £ t ]str[A 4 £<^£ t ] , (46) 
where we have defined the matrices 

= SisSji, (A 4 ) i;7 - = 6 a 6j2, (47) 

Aij = SuSj! - 6 i2 5 j2l B = diag(l, 2, -3). (48) 

Since the kaon has isospin / = 1/2, the AI = 3/2, K — > process vanishes. The 
amplitude for the (27,1), AI = 3/2, K — > rr matrix element is 



( 7r +| C )(27,l)(3/2)| ir +^ 



4«27 1 

.—m x m xz + — 



16 (-c 27 ,2 + 4c 27 ,24) m x m 2 xz 



+ 8 (c 27 ,4 - c 27 ,2o) m x m xz [m x + m 2 xz ) 



+ 8c 2 7 )7 m x m xz (2m 2 D + m|) 

where the logarithmic terms are given in Appendix [Fl 

For degenerate valence quarks, this amplitude simplifies to 



,7T 



| (27,l)(3/2)| jR: +^ )(4g) 



y\o {27 > 1){3/2) \K + ) 



dog.val. 



4a 



27 2 
2~ m X 



f 

8a 27 



1 + SZ X + 



5m 



x 



m 



x 



t" l~jt m x\ M(m x ) + 8£(m xd ) + U{m 2 xs ) - 3m x £(m£ 



+ 



P 



16 ( — C27,2 + C27,4 — C 2 j^ 2 q + 4C27 ,24) 



+ 8c2 7i7 m^ (2ra^ + m 2 s ) 



(50) 



VII. A/ = 1/2 WEAK MATRIX ELEMENTS FOR (8,1) AND (8,1)+(27,1) OPER- 
ATORS 

In this section we present results for AI = 1/2 amplitudes, which include K — > and 
K — > tt matrix elements, for operators that transform under the (8, 1) + (27, 1) represen- 
tation, and for those that transform under the pure (8,1) representation. We perform the 
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power subtraction explicitly through NLO in the chiral expansion and present the subtracted 
amplitudes. 



A. Definition of the C( 81 ) operator 

As shown in Ref. p|, there are two leading order [O (p 2 )} operators in the chiral symmetry 
group (8,1), with coefficients ot\ and a 2 . There are 13 NLO [O (p 4 )] operators relevant for 
this work. Note that there is an extra operator which only appears in the partially quenched 



case. In the full theory, by convention 



361 ] . operator 14 is absorbed into operators 10, 11, 



12, and 13 via the Cayley-Hamilton theorem. This is not possible in the partially quenched 
theory [ll|, Q. The (8,1) operators are 



=«istr [Ae^E^E+J + a 2 25 str [A 6 (A4 + E t + EM)] 
+ cgi^str [X 6 S 2 ] + c 8 i i2 str [A 6 S] str [S] 
+ c 8 i i3 str [X 6 P 2 ] + c 8M str [X G P] str [P] 
+ c 8 i i5 str [A 6 [P, S]] + c 8 i 5 i str [A 6 {S, L 2 }] 
+ c 8 i,iistr [A 6 L M SX M ] + c 8 i,i 2 str [A 6 L M ] str S}} 
+ CsL^str [X 6 S] str [L 2 ] + csi^str [X 6 L 2 ] str [S] 
+ c 81 , 15 str [As [L 2 , P]] + c 81i35 str [A 6 {L M , d u W^}} 
+ c 81 , 39 str [XeW^W^] (51) 

The (27,1) operators relevant for this section are given in Eq. (|4"U|) . 



B. K —¥ amplitudes 

The NLO expression for the unsubtracted K — > amplitude in the pure (8,1) case is 

(0\O^\K°) = ^a 2 (ml z -m 2 x ) + j(m 2 xz -m 2 x ) [2 (c 81>1 - c 81 , 5 ) m 2 xz 

+ c 81 , 2 (2^ + -l)]+<0|O (8 ' 1) |i,°> logs , (52) 

where the logarithmic terms are given in Appendix [G] For the (8,1) + (27,1) case, we have 



JO |0(8,D+(27,1)(1/2)| ^ _ |0(27,1) I ^ 



logs 



+ Y c ^H z - m x) 2 + (0\O {8 ' 1] \K°), (53) 
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where 



(0\OW\K°) 



logs 



«27 



2m 2 x -R x (m r ,) + m 2 x R n (m x ,m x ) 



+ 2m x R x (m rj ,m z ) 



{m 2 x ) + 



2m 2 z - R z (m rj ) 



+ m 2 z R r] (m z , m z ) + 2m 2 z R z (m v , m x ) 



£(m 2 z ) 



- m 2 R v (m x ,m x ] 



-m 2 v R v (m z ,m z ) + 2m 2 v R v (m x ,m z ) 
£(m 2 ) + m 2 x R x (m v )I(m 2 x ) 



+ m 2 z R z (m ri )i(m 2 z ] 



(54) 



Following the procedure given in Ref. 28J, we perform the subtraction of the power di- 
vergence in (8,1) and (8,1) + (27,1) amplitudes. In order to do this we require the amplitudes 
for e( 3 ' 5 > through NLO, given in Section HVl The ratio of (8,1) and (3,3), K — > matrix 
elements to NLO is 

( O |0(8>i)|K°) a2 



;o|e( 3 >3)|if°) 



0,(3,3) 

4 

+ 



B (m z - m x 



f 



a 



(3,3) 



m 



2ia^ 

x) [2 - 4i,b) m lz + 4l,2 { 2m D + m s)] > (55) 



where the transformed coefficients c' 81 1 , c' 812 , c' 81 5 [defined in Tabled] are linear combinations 
of the original LEC's, Csi,i, etc., and the Gasser-Leutwyler coefficients originating from 0(p 2 ) 
terms in the amplitude (0 |0( 3 > 3 )| K®}. The first term on the right-hand side of Eq. (155]) 
contains the power divergence, which is proportional to m z — m x to all orders in the chiral 
expansion. The remaining terms are finite, including the rotated LEC's c' 81i . Since the 
rotated LEC's contain a term proportional to a?2, it follows that the unrotated Cgi/s must 
also contain power divergences 37J. This was implicitly assumed in the work of Ref. 11|. 
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Table I: The transformation of the (8,1) LEC's (denoted in the text with a prime) under the vacuum 
subtraction process. 



LEC 


Transformed LEC 


c 81 1 


c 8 i i - (4a 2 /f 2 )(2L 8 + H 2 ) 


C81,2 


c 8 i,2 - (16a 2 // 2 )L 6 


C81,3 


C8i,3 + (4a 2 // 2 )(-2L 8 + ^2) 


C 8 1,5 


Cbi,5 - (4a 2 // 2 )^2 


C81,10 


csi^o - (4a 2 // 2 )L 5 


C81,13 


cgi.13 - {^a 2 /f 2 )L A 


C81.15 


CBi.15 + (4a 2 // 2 )L 5 



A similar expression exists for the ratio involving the (8,1) + (27,1) amplitude, 
'n |o(8,i)+(27,i)(i/2)| r^o\ f 

fr|e(M>|*») } ^* " + 21^) <oi^ 

+ 2 7i) (0|o(27 ' 1)|A " )l - 

+ ~ [ 2 ( C 81,l - C 81,5) m L + C 81,2 ( 2m D + m l) 

+6c 27j i(m 2 z - m 2 ^)] . (56) 



C. K — > it amplitudes with 2+1 valence flavors 



For the K — > 7r amplitudes, we first present the matrix elements of the unsubtracted 
operators, 



y\o^\K + ) 



T 2 

+ 



4a 



■m xz m x 



f 



2 2 
2 



.f 2 



-2c 8 i,imL - C 8 l,2^L ( 2m D + m s) ~ 2c 8i,3m 2 x m lz 



+ 2c 8 i i5 m 2 ^ [m 2 xz - m 2 x ) + 2c 8 i i i m x m^ + c 81>n m x m X g 
+ C8i,uwixm M (2m^ + m 2 s ) - 2c 81i35 m x m xz (m x + m 2 xz ) 



+ 8c 8 i !39 m^m 2 



(57) 
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<7T + |0( 8 ' 1 )+( 27 ' 1 )( 1 /2)| X + ) 



4d! 27 1 

T^^lrzlTT'X H o 

£2 xz f2 



A8c 27A m 2 xz ( 



m 



x 



> 



1 xz 



+ 16 (-c 27 , 2 + 4c 27i 24) m^m; 
+8 (c 27i4 - c 27 ,2o) ni x m xz (m x + m*J 
+ 2Ac 27i6 m x m xz (m x - m 2 xz ) + 8c 27 jm x m xz (2m 2 D + m 2 s ) 

■ |0(27.1)(1/2)| + ^+ | (8,1)| ^ (5 



,7T 



where the logarithmic contributions to the (8,1) amplitude are given in Appendix [Gj and 
the logarithmic contributions to the (27,1) amplitude are given in Appendix [Fl 

Given the coefficient of the power divergent term from K — > 0, it is possible to carry 
out the operator subtraction in K — > ir numerically. By CPS symmetry, the following 
subtraction removes the divergence to all orders in the chiral expansion, 



7T 



0(8,1; 

sub 



K+ 



,7T 



\° m \ K+ )- 2 ^ B o(m z + m 



7T 



0(3,3) 



(59) 



To NLO this expression becomes 

4a i 



7T 



^sub 



K 



f2 -m xz m x + — 



2 4i,i m L - 4i,2 m L { 2m D + m s) 
-2c' 813 m 2 x m 2 xz + 2c' 815 m 2 xz (m 2 xz - m 2 x ) + 2c 8110 m x m 3 xz 
+ c 8 i,um x m xz + c & i M m x m xz (2m 2 D + m 2 s ) 



- 2c 81t35 m x m xz {m x + m 2 xz ) + 8c 81 ^ 9 m 2 x 



y\ow\K+) 



(i) 

logs ' 



(60) 



where again, the d 81 i are the linear combinations of LEC's given in Table E Note that the 
subtraction eliminates the term in Eq. ( 1571) proportional to a 2 . The NLO chiral logarithms 
proportional to a 2 are also eliminated. The remaining logarithms are contained in the term 
given in Appendix [Gj this term is proportional to a±. 
A similar subtraction can be performed for the (8,1) + (27,1) case, 



7T 



o. 



(8,l)+(27,l)(l/2) 
sub 



K~ 



7T 



| C ,(8,1)+(27,1)(1/2)| 



a 



«2 



B (m z + m x ) In 



0(3,3) 



K~ 



(61) 
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Again, this subtraction removes the power divergences to all orders in the chiral expansion. 
To NLO the subtracted operator gives the matrix element, 



7T 



o 



(8,l)+(27,l)(l/2) 
sub 



A H 



4:0127 1 

■—m xz m x + — 



48c 27 ,i^L { m l 



m 



x 



) 



,2 2 



+ 16 (-c 27j2 + 4c 27i 24) m x m: 
+8 (c 27 ,4 - c 27j20 ) m x m xz (m x + m^) 
+ 24c 27i6 m x m :C2 (m^ - m^) + 8c 27 jm x m xz (2m 2 D + m|) 



+ <7r+ | ( 27 ,D(l/2)| K+)iogs+ ^- 



o 



(8,1) 
sub 



(62) 



For degenerate valence quarks, Eq. (|60|) reduces to 



7T 



o 



(8.1) 
sub 



deg.val. 



4 o 



(~ 2cg 1;1 — 2c' 81 3 + 2cg 110 + Cg^n 



- 4c 8 i i35 + 8c 81i39 )m 2 x + (-c 812 + csi,^) (2m|, + m|) 

+ < 7 r + |o( 8 ' i )|A' + >" o ; g ; aL ' (1) , 



(63) 



where the logarithmic terms are given in Appendix [G] For the (8, 1) + (27, 1), K — > ir matrix 
element, we have for the degenerate valence case, 



7T 



W8,l)+(27,l)(l/2) 



sub 



\ deg.val. 

A'+) 



4a 



/ 



2 7 2 



1 + 8Z X 



Srrix 1 



711 



X 



3 / 



4 m x \6£(m x ) + 8e(m 2 xd ) 



+ M{m 2 xa ) - 3m 2 J(m x ] 



1 

+ J2 



16 (-C27, 2 + c 27 , 4 - c 27j20 + 4c 27i24 ) m 



X 



+ 8c 27)7 m|- (2m^ + m|) 

deg.val. 



7T 



O 



(8,1) 
sub 



K~ 



(64) 



Note that for degenerate valence quarks the (27,1), AI = 1/2 amplitude is the same as the 
(27, 1), AI = 3/2 amplitude, Eq. ([50]). 
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VIII. FINITE VOLUME CORRECTIONS 



Incorporating the leading corrections coming from the finite volume used in lattice sim- 
ulations for the above expressions is straightforward. Here we assume that the time extent 
used to extract the above matrix elements is infinite, and that the only corrections come 
from the finite spatial volume. There are two classes of one-loop integrals that must be 
replaced by their finite volume counterparts. The first is defined in Eq. (IBlj) . and its as- 



sociated double-pole counterparts are defined in Eqs. flB2j) and (1B3I 
Eq. (IB 1 h by derivatives with respect to m 2 ]. As discussed in Refs. 
effects can be accounted for by making the replacements 



[these are related to 



3s 



391 ]. finite volume 



i(m 2 ) -> i(m 2 ) + -^—m 2 8 1 (mL) , (65) 

167T 2 

£( m 2) _ £ (m 2 } + J_<j 3 ( mL ) , (66) 

167T Z 

|( m 2 ) ^|( m 2 ) + ^ Ml^l , (67) 
167r z mr 



with 



5 3 (mL) = --^L [ m %(mL)] = 2 ^ K (\n\mL) , (69) 

n^O 

5 5 (mL) = m 2 ^- 2 [8 3 (mL)} = - ^(H^K^mL) , (70) 

n^O 

with Kq, K\ the modified Bessel functions of imaginary argument. 

The second class of loop integrals are more complicated, and are defined in Eqs. flB4j) 
and (lB5p . For these, we recall the technique used to calculate the above finite volume 
corrections. We begin with the finite volume Euclidean space version of Eq. (1B4I) . and 
apply the Poisson Resummation Formula (as discussed in Refs. j^s), Q). This leads to the 
following replacements, 

P (q 2 , ml, m 2 ) -> (3 (q 2 , m 2 , m 2 ) + —8 p {qL, m x L, m 2 L) , (71) 





mL) 




n 


mL 



/3 (g 2 , m 2 , m 2 ,) -» (3 (q 2 , mf, mf) - ^ 2 5p(qL,m 1 L,m 2 L), (72) 
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and the corrections 



5(3(qL,miL,m 2 L) = ^ / dk '"""^" u — r ^ r ^ ; 7^ — N21 , (73) 



DC 



fc sin(/c|n|) c^i + u; 2 



|n| co>iu; 2 [(gL) 2 + (oui + u; 2 ) 2 ] 
fcsin(fc|n|) (qL) 2 LU 2 + (2a>i + cj 2 )(^i + ^) 2 



0s(qL,m 1 L,m 2 L) = (m^) > / dfc n —3 — ■ r-y— , 74 



n^O 

where we have defined 



and where the function in Eq. (|74p is obtained by taking the partial derivative of Eq. (I73"|) 
with respect to m 2 . 

These formulas can be simplified as in Ref. 40|, but only in special cases (such as degen- 
erate masses). For the general case, it is more difficult to find an approximate expression for 
these finite volume corrections. 2 However, it is relatively simple to evaluate these expres- 
sions numerically at a finite number of points. Given a set of lattice data at a number of 
quark masses and lattice volumes, it would be straightforward to tabulate the appropriate 
finite volume corrections from the above formulas. 



IX. CONCLUSIONS 

This paper presents the calculation of K — > and K — > 7r amplitudes to NLO in PQ%PT 
with 2+1 flavors of non-degenerate sea quarks. Results are presented for both the AI =1/2 
and 3/2 channels, for chiral operators corresponding to current-current, gluonic penguin, and 
electroweak penguin 4-quark operators. The chiral operators are conveniently grouped by 
their chiral transformation properties; this work computes matrix elements of (8,8), (27,1), 
(8,1), and (8,1) + (27,1) chiral operators. The power divergent operator subtraction is per- 
formed explicitly through NLO in the chiral expansion for AI =1/2 matrix elements. We 
have also shown how to include finite volume effects through one-loop for the quantities 



considered in this work. These results are useful 
rently available 2+1 flavor lattice QCD results 



or studying the chiral behavior of cur- 



41], from which the low energy constants 



of the chiral effective theory can be determined. The low energy constants of these matrix 



2 One cannot apply the expansion in Ref. 40|, for example, because these integrals have three different 
relevant scales, as given by q 2 , m 2 , and m|. 
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elements are necessary for an understanding of the AI = 1/2 rule and for calculations of 
e'/e using current lattice QCD simulations. Electroweak penguin K — > txtx matrix elements 
can be constructed to NLO in xPT using the formulas presented in this work, allowing the 
convergence of the chiral expansion to be studied. This will serve as a useful cross-check for 
other non-xPT methods such as those proposed in Refs. j|3, 
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Appendix A: ISOSPIN DECOMPOSITION 



The operator that governs the transition K — > ir can have either isospin 1/2 or isospin 



3/2. In typical lattice calculations, these two processes are calculated independently |28l. |43| . 
since the operator with isospin 1/2 mixes with a divergent lower dimensional operator, which 
must be subtracted. The isospin 3/2 amplitude does not have this complication. Therefore, 
we calculate the amplitudes for these two processes separately, making use of the amplitudes 



for K + — > ir + and K { 



7i 



We define Ai + = (ir + \Oi\ K + ), where Oi represents some AS = 1 operator with both 
isospin 1/2 and 3/2 components, and Aio = (tt° \Oi\ K°). If we decompose the operator Oi 



by isospin, Oi = O. 



(3/2) 



o 



(1/2) 



then we have for the matrix elements, 



M + 
Mo 



7T 



7T 



o 



(3/2) 



o 



(3/2) 



K~ 
K° 



71 



^(1/2) 
i 

0(1/2) 



K 
K° 
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Given the relevant Clebsch-Gordon coefficients, 



we obtain the result, 



1 

(n + 


^(3/2) 


\ 

K+j 




£,(3/2) 






0(1/2) 


K+\ 


(n° 


0(1/2) 





71 



71 



o 



(3/2) 



o 



(1/2) 



K~ 
K~ 



V2 
2 ' 



-V2, 



- 3 [M + 



2M. 



y/2M 
-V2M 



(Al) 



(A2) 



(A3) 



Appendix B: LOOP FUNCTIONS AND RESIDUES 

The following loop functions are used throughout this work, and they are regulated using 
dimensional regularization in the modified MS scheme. For single-pole mesonic loops, we 
need 

d d p 



£(m 2 ) 



[cf. f 2 A (m 2 ) in Ref. 



10 



lim 



1 2, / 
m In 
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111]). We also need the double pole expression 
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where the minus sign is chosen to be consistent with the form of Euclidean i (m 2 ) in Refs. 39, 
44l| . 3 Further, we will sometimes need 
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The two-point loop-function encountered in loops with strong-weak vertices and only a 



Note, however, that our definitions of I and I differ from Refs. 
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44| by a factor of 1/16tt 2 
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single pole is denned as: 



(3 (q 2 ,ml,ml) 



d d p 



(2%) d (p 2 — m\) ((p + qf — mfj 

— 2 / dx {l + In [— x(l — x)q 2 + (1 — x)m\ + xm\?\ — In (/x 2 ) } 
7TJ Jo 



(4tt) 



(B4) 



Note that we always have q 2 = (m xz — mxf for K — > 7r amplitudes. This function is 
proportional to the B function defined in Eq (A2) of 10J. Similar loops with double-poles 
require 
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To simplify the expressions, we use the notation for the residues arising from disconnected 
meson propagators, 

(ml - m 2 D ) (m 2 x - m 2 s ) 
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Appendix C: ONE-LOOP WAVEFUNCTION AND MASS RENORMALIZA- 
TIONS 



The necessary wavefunction renormalizations needed for the one-loop amplitudes are 
given in the 2+1 flavor case by 

SZ XZ = -2^ + - (^f) , (CI) 

J 6 V / /logs 

where we have separated the terms in this way because the first term on the right-hand side 
of each of these equations is the one- loop correction to the bare decay constant / appearing in 
the tree-level expression for a given weak matrix element. This first term contains both NLO 
logarithmic corrections and Gasser-Leutwyler constants. It may be useful in chiral fits to 
lattice data to absorb this correction to the decay constant into the tree-level expression, and 
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the above formulas make this convenient. The second term is proportional to the logarithmic 
corrections to the decay constant alone, without the Gasser-Leutwyler constants, 
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For the degenerate mass case, these reduce to 
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Additionally, we need the one-loop corrections to the meson masses squared jsj], 
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Appendix D: LOGARITHMIC CONTRIBUTION TO (3, 3) K -> ir MATRIX EL- 
EMENTS 



The logarithmic contribution to the (3,3), K — > ir matrix element for the 2+1 non- 
degenerate case is, 
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Appendix E: LOGARITHMIC CONTRIBUTION TO (8,8) K -> ir MATRIX EL- 
EMENTS 



The logarithmic contribution to the (8,8), AI = 3/2, K — > ir matrix element in the 2+1 
non-degenerate case is 
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The logarithmic contribution to the (8,8), AI = 1/2, K — ► ir matrix element in the 2+1 
non-degenerate case is 
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For completeness we include the chiral corrections to K — > txtx at physical kinematics for 
the electro-weak penguin operators. In the full theory at physical kinematics, the logarithmic 
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contribution to the (8,8), AI = 3/2, K — > rnr amplitude is 
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and the logarithmic contribution to the (8,8), AI = 1/2, K — > tttt amplitude is 
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Appendix F: LOGARITHMIC CONTRIBUTION TO (27,1) K -> vr MATRIX EL- 
EMENTS 



The logarithmic contribution to the (27,1), AI = 3/2, K — > n matrix element in the 2+1 
non-degenerate case is 
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The logarithmic contribution to the (27,1), AI = 1/2, K — > ir matrix element in the 2+1 
non-degenerate case is 
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Appendix G: LOGARITHMIC CONTRIBUTION TO (8,1) K -> AND K -> vr 
MATRIX ELEMENTS 

For the (8,1) case, we separate the logarithm terms which are proportional to a\ from 
those proportional to a 2 - 
For K — > 0, we have 
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m x 



+ 



2m| - Am 2 xz 



m xz - m x 
R x {m v )I{m 2 x ) + 
3mx{rn 2 z — m 2 x ) 



3mx 



m xz - m x 



P(q 2 ,ml,m 
Rz(jn^)l{m 



2 ) 

xz) 



4; 



4m, 



Rzim^piq 2 \m 2 z ,m 2 xz ) \. 



m xz - m x 

These formulas are simplified enormously in the degenerate valence case, 



/„+ 10(8,1)1 ^+\de g .val.,(l) 
\ / logs 



5Z 



x 



Aa i 
7 



Aa>i 2 I 



5m x 



rn 



x 



+ [m 2 +m 2 x ] (3(0, m 2 , m 2 x ) 



R v (m x ,m x ) 



2 - R v (m x , m x ) 



l{m\) - 10£(mL) 



— bi(m 2 xs ) — 2m 2 x Rx(m v )£(m 



+ 4 



-m\ - Rx(m v ) + m 2 x R v (mx, m x ) 



Arn 2 x) , 
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U + \ o^\K + ) dcg '™ l " {2) 

\ I I / logs 



8a 2 2 J 2 
3f 4 x \ x 



R x (m v )i(m 2 x ) - 2R n (m x , m x )/3(0, m 2 , m 2 x ) 



-1 + R v (m x ,m x ] 
R v (m x ,m x )£(m 2 ) 



i{m 2 x ) + 2£(m 2 xd ) + £{m 2 xs ) 



+ 



2m 2 x 1 - R rj (m x , m x ) + 



(G8) 



Appendix H: ERRATUM 



We note here some corrections to the works of Refs. 



10| and Al 



energy constants for the (27,1) operators have the wrong sign in both Ref. 



The values for 7, appearing in Table I of Ref. 



the NLO low 



10| and Ref. (ll| 



. 3 



10] (and again in Table I of Ref. 11 j) should 



(and again in Eq. (16) of Ref. Ill]), the 



have the opposite sign. In Eq. (16) of Ref. 

(8 1) (S 1) 

operators 0§ ' and 0{§ should have opposite sign to be consistent with the signs of the 
LEC's e r 5 and e r lb appearing in the amplitudes presented in those works. We make these 
corrections in the current work. Since the LEC's are not known, an incorrect, but consistent 
normalization of them (including an incorrect sign) does not alter the procedure of using 



the formulas of Ref. 



llj to construct K — » tttt from K — > 7r and K — » matrix elements. 



Therefore, these corrections make no difference to the conclusions of these works that it is 
possible to obtain all of the LEC's needed to construct K — > tttt matrix elements through 
NLO in %PT from lattice calculations. 

Additionally in Ref. 11], there is a typo in the (8,8) K — * tttt matrix element formulas. 



The coefficient of the second term in Eq. (42) should be j^jz, and the coefficient of the second 



term in Eq. (43) should be 



12/ 



Sufi' 



The corrected versions of these equations are given in 



the current work as Eqs. (1381) and ( 1391) . These corrections also do not alter the conclusions 
of Ref. ll|, but are necessary to construct the correct K — > tttt matrix elements. 
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